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ABSTRACT 



The solutions of the phase equation for a potential with the 
asymptotic properties of an atomic polarization potential are 
studied with the intention of developing a computer program for 
direct integration of the phase equation for the s-wave phase 
shift of an arbitrary potential. Such a program is developed and 
its limitations are discussed. In addition, an appendix is devoted 
to the variational principles of Kohn and Schwinger for calculating 
phase shifts. 
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THE 


K 


THE 


n Q (r,K) 


THE 


n o (K) 


THE 


u(r) 


THE 


U(r) or V(r) 


THE 



RADIAL COORDINATE IN BOHR RADII 
ENERGY IN RYDBERGS 
PHASE-FUNCTION IN RADIANS 
PHASE SHIFT, n (»,K) , IN RADIANS 
SCHRQDINGER RADIAL WAVE FUNCTION 
POTENTIAL FUNCTION IN RYDBERGS 
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I. INTRODUCTION 



The Quantum Theory of scattering processes is of basic import- 
ance in studies of atonic, nuclear, and elementary particle physics. 
Analysis of Atonic collision processes is required in such diverse 
applications as gas dynamics , chemical reaction theory, and stellar 
dynamics. Nuclear Physics relies on scattering of various particles 
fran nuclei for a great deal of the quantitative information about 
nuclear processes while in Elementary Particle physics, very little 
other information is available except for the various scattering 
events that are observed. 

There are many approaches to the scattering problon, but one of 
the most useful (especially in low energy applications) is the method 
of partial wave analysis. The method was first introduced by Faxen 
and Holtzmark based on a method developed by Lord Rayleigh. Rather 
then delve into the basis of the theory, it will be assumed that the 
reader is familiar with the treatment of Schiff [Ref. 1] or Merzbacher 
[Ref. 2]. 

A particularly interesting procedure for calculating scattering 
phase shifts is the "Phase Method." Originally introduced by Courant 
and Hilbert [Ref. 3] and applied to scattering by Morse and Allis 
[Ref. 4] , the method has recently been revived and expanded in a 
series of papers by Calogero [Ref. 5, 6] and his new book [Ref. 7] . 

This paper will investigate the application of the phase equation 

to a potential which has the asymptotic properties of an atomic 

4 

polarization potential (1/r ) . The method employed is a direct 
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integration of phase equation (11) . In the appendix, the variational 
principles of Kohn and Schwinger are treated and are suggested as 
alternate methods for computation of phase shifts. 
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II.. THEORY 



A. THE PHASE FUNCTION METHOD 



The phase function method is a method of reducing the second order 
linear Schrodinger Equation to a first order, non-linear Riccati equa- 
tion. For s-wave scattering, the solution for phase shifts reduces 
to the solution of the differential equation 

(1) rj'(r) = - [ V(V)/k ] t rj) 

To show hew this equation canes about, it is convenient to approach 
the problem frcm the standpoint of A. Ronveaux [Ref. 8] . He starts by 
considering a second-order homogeneous differential equation with two 
perturbative terms , (X) and (X) ; 



{2) U." ■+ fKx) u' + Q(x)u =■ \J t (x)u. f V z (x) u' 

with u(0) = \i and u' (0) = v. This can be written in the form of the 
following system 

(3) X'(X) - A(*>X(*> + VCx)X(x) 

where / K , / \ 

X ( x ] - M X(o) -- 

A(x) -I O i \ VO) -to O \ 

(-Q -P/ ( V, vj 

The unperturbed solution, X q is, 

( 4 ) XJx) = W(X) C 
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where C 



is a constant vector, and 



U, 



W(x) - , , 

ix t 



lA^ 

U- 



and u^ and are two linearly independent solutions. The solution of 
the complete system may then be written 



(5) 



X(x) = W(x)CCx) 



Using the Lagrange Method of Variation of Constants, then 

C'(X) ■ 1 W«)!'Vw W(x) C(/) 



or 



where 



C'(x) -- MO)C60 
C (o ) - ^ Wfo)| 0* j 



If we now define the phase function 



(6) $(x) - C 2 C *)/C,(<) 

and differentiate, we find that 

(7> s'(x) = n 2s M + (x) - a),, wj so) - m, 2 o)sVx) 

Evaluating M gives the phase equation 

= U^U)] + [u 2 a)SOr) 

( 8 ) 

X I v;w[u,W t Kjdf) SW] + V 2 W[u,M i-w/MSWj 

where det(W)^ 0 and 




S(o) = 



1! 


1 " 




[ u,(o) y] 


1 


[J 


U z ^o) y juj 


f * i 


f y 


i 
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Given the general form one now approaches the problem of s-wave 
phase shifts. The Schrodinger Equation (in natural units, h " C - I ) 
is 

(9) Ul'X V) -f K l u(r) - V(r) tur) 

choose 

U r yA. r. O 

l)J io\ = V 

Then the general solution is 

U(r) ' C, S/n fcr - C 2 coi kr 

^ C t (r ) [s in Kr - $(r) Cobkr] 
or 

UCr) ~ C, (r) coi r ? 0 ( r, K) S'* ( Kr t ty) 

where 

s(r)= C 2 (c) /c, (r) = - -fan y^Cc^ k) 

Then, letting V, - V 2 = 0; U,- Smkr, u z - cos/< r, dt'f U/; K the general 

phase equation becanes 

(10) (d /cl r ) "fan = - [ V(r) //cj S'n Icr -f fa n co $> /c'r j 

with "fan (o,K)z o. Noting that (d/dOfon^ : S<fC * y ]o( d *] 0 /d f ) 

it beccmes apparent that 

I 11 ) (d/dr^o = - [ V(r)/K j si-i 2 (l<r + ) 

B. THE POTENTIAL FUNCTION 

The object of this paper is to investigate the properties of the 
numerical solutions for the Phase Equation (11) when applied to a 
potential which has the asymptotic properties at infinity of an 
Atonic Polarization potential, i.e., 

U(r) > (const)* 1 /r ^ 

r — ■? 
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In order to simplify the analysis, an additional condition was imposed 
on the potential in the form of a requirement that as r-> O , the 
potential should approach zero as f , i.e., 

U(r) — - — — —> (consf)xr 
r •— *•> o 

The function chosen to satisfy these conditions was 

U(r) -- 1 fat ^e' + ? jU 

] r—?oo 

Jo 

The asymptotic properties of this form are obvious, however, it is 
difficult to calculate with. In order to simplify the problem further, 
the substitution T - X r was chosen which leaves 

(12) UCr) ~~ ridyx^e rX 

Jo 

It is no longer readily apparent that U (r) has the correct asymptotic 
properties, but it may be easily verified by integration and power series 
expansion. The success of this form is shown in Fig. 1. 
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V (r) 




Figure 1 

Plot of the Potential Function 
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III. CALCULATIONS 



Having derived the phase equation (11) , one now wishes to devise 
a method to put it to practical use. Calogero develops a variational 
scheme for the phase equation [Ref. 9] while other authors deal with 
the equation in other ways [Ref . 10 , 11] . An excellent review has 
been prepared by Babikov [Ref. 12] which provides a good general 
introduction to actual calculations and an equally good bibliography. 
Owing to the fact that the equation is non-linear, the various ap- 
proaches have dealt with sane form of approximate solution. This 
paper will deal with the relatively simple problem of constructing 
a program to directly integrate the phase equation using numerical 
approximation only. 

The particular method chosen was a fourth order Runge-Kutta 
Integration scheme to be implemented on the IBM-360 time sharing 
system. General discussions of the numerical method may be found 
in numerous places [Ref. 13, 14], and so only the necessary portion 
of that theory is presented here. 

The Runge-Kutta method proposes to solve the equation 



subject to the initial conditions X = X when Y = Y . The iterative 
J o o 

solution process proceeds via the following equation? 



(13) 




(14) 




(o 



where 
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Km r ^ ( Xrj * ^ f 1C3) 



and where / r)+| = f )i . The progr arming needed to provide a 
solution of the equation (11) which gives the exact phase shift, 

( t<) - k) / then reduces to programming a simple iterative 

procedure. The program developed is printed at the back and is in 
such a form that it may be run either on the time-sharing system for 
quick approximations or directly on the computer system in the regular 
batch processing. A second version of the program used to obtain a 
print out of the values of the phase function at various r values is 
also included. The use of the second form will became obvious as 
the discussion unfolds. 

It now becomes necessary to consider the actual running of the 
program and its particular limitations. In order to properly investi- 
gate the program, the Yukawa potential, 



is convenient to consider since Calogero treats the case with strength 
parameter, g = -10, m = 1. The output produced by the second program 
is graphed in Figs. 2 and 3. 

The most obvious feature of the graphs is the presence of "levels" 
which arise from the presence of possible bound states in the potential. 
Calogero discusses the connection of this with Levinson's theorem and 
shews that this structure in the phase function can be used to locate 
bound states. This very feature though poses a distinct problem to 
the evaluation of the exact phase shift. Owing to the presence of this 
structure, a simple test for the end of the integration that depends on 
the phase function tending to a constant value at infinity is useless. 
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Figure 2 



Plot of Phase Function for K = 0.1, Yukawa Potential 
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Figure 3 



Plot of Phase Function for K = 1.0, Yukawa Potential 
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For this reason, the integration must be carried out past any possible 
structure before a cutoff of the above type is allowed. In order that 
the program will do just that, an integer counter is established and 
required to be larger than a certain lower limit which is designed 
to carry the integration over the structure. The procedure for doing 
this is as follows; the second form of the program is first run with 
the intention of finding the structure and determining the value of DR 
which gives the best results. Having found the structure, the counter 
may be set appropriately. By use of this artifice, the program is 
enabled to reproduce the results given in Calogero's book to within 
a few percent. It should be noted that another way of treating this 
problem is to establish a lower limit on the basis of the value of r 
reached. This approach would be superior if there were a method for 
predicting the position of the last flat region in n Q . The counter 
method was chosen only on the basis of a desire to know the number 
of iterations used, for computer timing purposes. 

A further pecularity with the direct integration routine is the 
sensitivity to the x increment, DR. In investigations with different 



* * a * 



40 



values of DR, it was found that when 



i 



DR is chosen too large the structure 






t 



at small r tends to destroy the ac- 
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curacy of the routine by adding 



large terms over some of the flat 
regions in the phase function. An 




example is for K = 0.1 and dr = 0.1. 



Figure 4 



The result is n = 44.1 in 110 steps 
o 

as shewn in the diagram (Fig. 4) . 
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The true result (in agreement with Calogero's book, p. 22) is 
n Q = 6.24 in 38 steps found by taking dr = 0.01 or less. Now the 
utility of the second program for preliminary tests become apparent 
since it allows us to establish the number of bound states for the 
potential and the proper increment and the lower counter limit in 
just a few test runs which may be quickly run on the computer time 
sharing system. 

The program also has a very important limitation in that the 
accuracy decreases sharply for a small k due to numerical round off 
errors arising from the necessity of dividing by a small number. The 
situation may be improved by choosing a smaller increment, DR, but 
after a certain point this procedure tends to become self defeating 
because decreasing DR increases the program running time. A test run 
on a Yukawa potential with g = -2, K = 0.1, and dr = 0.1 gives a close 
answer to the zero energy phase shift, but if the value of k is reduced 
to K = 10 3, the accuracy in the calculated value of the phase shift 
actually decreases. 



21 



IV. RESULTS 



The final problem is the application of the programs to the 
polarization potential discussed earlier, eqn. (12) . Figures 5 
and 6 detail the results of the explicit calculation with the 
second program. They indicate that a value of DR of 0.1 will carry 
the program nicely through to the asymptotic region if the lower 
counter limit is set at 200. The lack of structure in the graphs 
tends to make this requirement useless though and it becomes pos- 
sible to set the limit at 10 iterations just to let the program get 
started before the asymptotic cutoff procedure is allowed to function. 
The lack of structure just carmen ted on tends to indicate a lack of 
bound states in the potential as might be suspected frcm the very 
weak nature of the function. Tests could be done by imposing a 
strength parameter on the potential, i.e., 



which could then be adjusted as suited the investigator. By means 
of such an artifice, bound states might begin to appear when the peak 
of the potential is high enough. In addition, increasing the strength 
of the potential would increase the magnitude of the phase shift. One 
might expect that such structure would appear for values of r less than 
1.0 and that the way to locate such structure would be to decrease the 
value of the increment DR. Though interesting, such an investigation 
does not have any apparent physical use since the only physical feature 
of the potential is the asymptotic dependence at large r. Thus further 
discussion of this point has been emitted here. 
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Figure 5 



Plot of Phase Function for K = 0.1, Polarization Potential 
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Figure 6 



Plot of Phase Function for K = 1.0, Polarization Potential 
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The most interesting feature of the calculation frcm the stand- 
point of possible application of the program to real, physical 
potentials is the necessity of carrying out the calculation to larger 
values of the radial coordinate than was the case for the Yukawa 

potential. The reason for this behavior is obvious, the polarization 

4 

potential only drops off as fast as 1/r at infinity. This leads to 
the observed difference in the convergence of the scheme and probably 
to the reason why the only places that the Phase method is applied in 
the literature are in papers dealing with nuclear type potentials which 
may be assumed to fall off very rapidly. To get some idea of the magni- 
tude of this effect, the graph in Fig. 2 may be examined in comparison 
to the graph in Fig. 5. The plot for the Yukawan potential has sta- 
bilized by r = 1 while with the polarization potential, the plot is 
just stabilizing at r = 25. These observations hold great import for 
any application of the method to atomic like potentials since the matter 
of time available for computation is often very critical. Even on the 
IBM-360 time sharing system, timing for the program had become critical 
when increments less than 0.1 were chosen. In addition, output disk 
space became a problem. The problem of timing could be expected to be 
considerable worse on smaller computing systems. 

Having discussed the phase function and its properties, use may 
nc m be made of program 1 to obtain a plot of n o vs. energy as in Fig. 7. 
The most apparent feature of Fig. 7 is the 0.78 radian peak at k = 0.36. 
It is a real peak as opposed to an apparent peak that could be produced 
by a poor choice of the increment at such a low value of k. An inter- 
pretation of this peak will not be attempted at this point but a com- 
ment is appropriate on the behavior of n Q (K) on the low energy portion 
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of the peak. It seems to drop off drastically, possibly to zero. 
Belcw k = 0.1, no data taking was attempted due the necessity of 
integrating too far with too small an increment. 
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Figure 7 



Plot of Phase Shift vs. K for Polarization Potential 
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V. CONCLUSIONS 



The program developed for the phase equation integration has been 
tried on a test case and applied to a second, semi -physical case with 
somewhat the expected results. This program has in these two cases 
performed well as long as the precautions mentioned were taken. It 
should be noted that no error analysis of the routine was attempted 
due to lack of time, however, the comparison of the results on the 
Yukawa potential for g = -10 with the results published by Calogero 
in his book [Ref. 7] , tends to point to an accuracy that is within 
a few percent. Whether this result is valid for any potential other 
than the Yukawa potential and in particular whether it is equally 
true for atonic level potentials is a question that remains unan- 
swered at this point. It is expected that simple numerical errors 
arising from the Runge-Kutta integration scheme and the ordinary 
computer errors of round off, etc. will decrease the accuracy of the 
calculation. 
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APPENDIX 



(A.l) The Kohn-Hulthen Method 

The most frequently used variational procedure is the Kohn- 
Hulthen method [Ref. 15, 16] which is based on the stationary property 
of the functional 1 where 



(Al) 



I 



= 2 rn f tp(e- H) 4 * dr 
J 

and is sane function of the phase shift determined frcm the normal- 
ization of 6 as r — y o<? . This paper will snploy tangent normalization 



throughout , i . e . , 

(A2) f) > (t/n) 2 + V} x C VT r ) 






where ^ is the phase shift for the partial wave and jj? and 

are the spherical Bessel and Neuman functions respectively. Consider- 
ing the special case of s-wave scattering, 5 = 0, then (A2) becomes 

4 7 2 [ Sinl<i-/|<r + <wKr/i rrj 

= (•&> 



cos 



*?<> 



5 X r t 
Kr 



(A3) 



61 



With this form for 'P , one new considers the explicit variation of (Al) 

j<56(v 2 + K z - u)^ ar t- ^(^ + k 2 - u)^<f r 

Since for s-wave, ^ has no angular dependence, 
so that the second integral in (A3) becomes 



! r 



-oo 



- ^Jl^meae y4oc*-o)^f 



29 



Integrating by parts, this gives 






oo 

dr • 




( l|) r ‘ dldf') 


i 


1 dr 


1 dr 


oo ( 






U r 


4> r 2 «! (<5<P) 


i i 


dr | 


u ^ 



- cL^ r 1 d ^ J*/') 

t|r 7 

+ ^d + C^-oW-'d' 



4 






oc 



- ^»r \c\r 4 - 3 ^r l cH 

ctf d«~ dr 



2^ {ok£"«H<l' 



4 w _„' \ o 

4s 



Thus (A3 ) beccmes 

SI = 2 ^‘K'v^hc*- U ) 4*cl 



co 



? + t-/ r. 



r z (^ cL£^~ d4 *cL^ 

( Jr dr - /. 

The last term, upon application of the boundary conditions, reduces 
to- ^GartVJo/O so that 

04 , £(I+ +•»>/*) ' 

But the Schrodinger equation is just 

( V Z 4 l< 2 - 0) P - O 

so that 

^ <^(I 4 4 O VI Y]o /fc*) - O 

The stationary functional, [A] = 4o«v^/je 7 j_ s given by 

< a6 > m -■ x + 1 5 



\ -t ( ij> ( E ” H) f ci 
\ 



r 
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The trial function 



ion, ^ , 



is chosen to have the form 



(A7) 



<p 1 X t (r) 

__ \ 

where and the Xj are sane convenient basis 



functions, that rapidly approach zero as r-? 0 # . The asymptotic part 



is chosen so as to shield the — singularity 



of the wave function , 

in cos kr/r as r-; 0. It is evident that has the correct asymptotic 
form as 1 



(p s (lJji ) * [si* kr/Kr - A 1C (*>•« kr/i< 2 r * - C oi kr /«r)J 

> ( *■/ 77 ) 2 [ ( s i n If r her) +A (coslfr/r)] 



r-? «= 



The Kohn-Hulthen principle then states that the variation of $ ] with 
respect to all arbitrary constants, i.e., 

(A8) £[\] £(\ + 1 ) * 0 

produces a value of [A ] that is second order accurate in 2 A • The 
explicit variation is done by writing the functional in the form; 



(A9) 


[A] - 


Iq 


C j ^2j 


Cj 2 j + 3 


where 


= 


2m 


^ X/ ^ - H ) - 






= 




(£• h)*P 






ft -- 


2 J23 

4? 


^ tf(e- m) <P 
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Variation of the n constants, then gives the system of linear equations 

(aio) 7 M ij Cj - - 

J 

Since the trial function is linear in \ , 

Hi ■- i + X fc/' 1 

and since the are independent of /\ 

+ A c; 

This leaves the equations: 



CO 



(All) 



i Hi, . ~sr 



I Mfi c/ ,N » - e. 



C/A 



Solving for the constants, 

= W c 4 i W| 4 /\ 4 4- \ B, + A Ba A 

where 

, <T CoA aj , 9 5- r , fo > n (0) 

c i Mij + 2 2 C i Rj 

* 

Co) - ^ ~ - - C/A.^, <■<*> 



\aJ, 



= 2 c“/7i, c,"> t 2% (.■/ 2; t 2 z c, /?,- 
jj ' J ' ' 



W- 



. pr^r * *z*r*r 

L 

Then a final variation with respect to \ gives, to first order 

( W| 4 B, 4 | ) 



A" 



2 C I - 1 S 2 ^ 

and the second order accurate, stationary functional is 



(A12) 



[A] = Vo i Bo - <^v, <- 6, i- n 

^ (w z + & L ) 
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(A. 2) The Schwinger Method 

A less frequently used but more elegant variational method is 
that due to Schwinger [Ref. 17, 18, 19, 20]. It gives the stationary 
value of the functional [Kco+v^] 

directly. The problem is to solve 

the radial Schrodinger Equation. 

(AI3) I-Vu, + \ K 2 - SLSjlO 1 u,(r) -■ U(f)u r (r) 



l 



dh 



To do so, one solves for the Green's Function Ci^Cr )f 0 ) given by. 



(A14) 



it C)£ ( r* I O + 



dr 3 



k 1 - Jg..C S L± J. 



Cij(rlr 0 )- - c )(r- n) 



which gives 

G ne(rlr 0 )~ ~Krr 0 (ifr ) (k r c ) r < r c 

(A15) * 

- - k r r c j A (. \< r 0 } y\ a C kr ) f a < f 

The integral equation for U f (r) is then 






(A16) 



Ojflr ) - Krjjttr) 



d r c C l0 ( r / r c ) UC r c ) u e ( rj 



For large r, the asymptotic forms are 

l< r jp ( \( r ) 7 c o < (hr - i ( + i)n ) 



k r ( \c r ) 



? S m (wr - ; i)u) 



therefore, the asymptotic form of U B (r ) may be written as 
Up( f ) Co t [ K r - £ <JHI)T,] 



(A17) 



. Oo 



-h [iCr - (S+Ot] r « r o I* < * r 0 ) (J ( r c Wa ) 

Jo J 

But, this is to match to the asymptotic form 

U e Cr ) — 7 coi [tr - ~ (e+i)" J -f si* fkr - -j (-P+ i)77_/-fa>o 
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and so one sets 



^ ^ - \cVr 0 r e U( r Cl ) u, < O 

J O 

Inserting the integral form for jtf(K') in terms of ty(r) frcm (A16) 



.or- 



4av'T| fi : i \ c) f 0(^VIO \ U? < * ) +■ ICiCMr*?) U(fo ) ^ C r t .) cl f c 



therefore , 



|C fa n t U(f ) u c 2 (f) -t \ d <0 OC r ) <Jg (■ f ) G\Cr Ko ) UCf^) o e (f c ) 



dividing new by 



-fa 






(A18) 



KCO^Ij - 



^ cW OC f ) Cr ) t jj”dr j d'e UC f ) lr ^0( r o^ t v<0 



i <J r f U ( r ) u p ( r ) i e ( t r ) 



For s-waves this becomes 



(A19) 



KctJrj 0 r 



, d ( U ( c ) o \ r ) - 






c\r 0 0(f) vie) (p(r\r 0 ) 0(r o ) v Cr c .) 



~ \dr UC f)i’(r)sin Irr 



J 



where G]( r lr 0 ) is new: 

C-](r I r 0 ) - s i n k r c o s k t~ 0 r < f u 

N 

r i. Sin l< coi It f 

K 

In order to show that (A19) is in fact stationary to second order for 
the exact wave function, consider 



(A20) 



Kcoi 0 - £ k, C c 4 ( c t <9 u 1C c o t f ^ 
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Expanding and collecting terms, we have 



£>{.Kcoi UCr)o(r)^oCr)- ^cir 0 (rfinr^) 0(^o) u( ' f c ) \ 

4 2 ^cir U(r)dv [uCr)~ 
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ccmbining the terms the numerator for the expression will be 
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- ulr)- ^Jr 0 C-jCr \r 0 )o(r c )v‘r tt )J jfejclr U(r)vO)s, 0 l<j 

+ aj^dr U(r)So«i sir>Hf J * ' jjdr 0 (r) u(r)$mkfjj -f (DCS 1 ) 



However 



jdrCXVWrj u <r)- \^i 0 C-)(rlr 0 ) 0(r o )o(r 0 k c o'f ^ jc\rCXr)v(r)smkrJ 



(A21) 

but frcm (A17) this is 



.a? 



dr U(r) o( r > 5 1 * /< k 



After factoring out jj| jcli 1 UfrluCr) sirtrj , the numerator beccmes 

[ic j^dr OCr)^(r) s<»kr - r kr 

4 2 jd r U(r) &v(r) £<->(r) - staler - ^d ^ £j <7 / £ ) U (r 0 ) u(r 0 )jj i ) 
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However , combining terms , this is , 



O(r)o(()swl , r | ^ 2. (J(r)£u(r) 

+ «6 l ) 



d r 0 0)1 r )r 0 ) U(fij>u<^) J 



But from (A16) , for u(r) the correct wave function, 

C* 



o(d - - Ur 0 G\(r\K)o(r 0 )vir o '> ~ O 

id z 

Thus, neglecting terms of order <5 / the numerator is zero and hence 

&[\i o 

Therefore, we see that (neglecting second order variations) the expres- 
sion is stationary when u is the correct solution of the Schrodinger 
Equation. 

New going back to equation (A19) and inverting 
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where 
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The problem remaining is to solve for and the C s . To 

do so, define 

(A26) Q: I c i c j ) £ o 

but Q is obviously stationary and so, 

(A27) r 2 - A///) = 0 

I 

In order to have a solution for the G$ , the determinant of the matrix 
of coefficients must be zero, i.e., 

(A28) Hij " AJij ) - O 



However, it is a simple matter to show that the matrix N is of rank 
one, i.e., there are no non-zero sub-determinants of the matrix N of 
order 2, and further, it can be shewn that such a system has only one 
nontrivial solution for £ K 4 £ *1 1 > 



(A29) 
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Thus one finds the solution for j l(i 0»t yj without recourse to a 
solution for the C5 . This result is shown in Kohn’s Paper [Ref. 15] . 
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